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Abstract 

■ We study the problem of optimizing nonlinear objective functions over bipartite matchings. 
O^l . While the problem is generally intractable, we provide several efficient algorithms for it, including 

a deterministic algorithm for maximizing convex objectives, approximative algorithms for norm 
U: minimization and maximization, and a randomized algorithm for optimizing arbitrary objectives. 

o" 

1 Introduction 

Let N := : 1 < i,j < n} be the set of edges of the complete bipartite graph K nn . In this article 

we consider the following broad generalization of the standard linear bipartite matching problem. 

o 

Nonlinear Bipartite Matching. Given positive integers d,n, integer weight functions w 1 ,...,w d 

■ on N, and an arbitrary function / : Mr — ► R, find a perfect matching M C N maximizing (or 

o ' 
\o 
o 



minimizing) the objective function f(w l (M), . . . ,w d (M)) where w k (M) := ^2{w k (i,j) : G M}. 



Identifying perfect matchings in K n>n with permutation matrices and weight functions with integer 
matrices in the usual way, the problem has the following nonlinear integer programming formulation: 
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>! max{f(w 1 x,...,w d x): ^Xjj = 1, ^Xjj = 1, x € N" xn } 

where w k x := Y17=i S?=i w ij x i,j f° r ^ = 1> • • • > and where N stands for the nonnegative integers. 

The problem can be interpreted as multiobjective bipartite matching: given d different linear objective 
functions w 1 , . . . , w d , the goal is to maximize (or minimize) their "balancing" given by f{w l x, . . . , w d x). 
The standard linear bipartite matching problem is the special case of d = 1 and / the identity on R. 
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Beyond the intrinsic interest in studying the above natural nonlinear extension of the standard 
matching problem, it is interesting to consider it in connection with its various variants and relatives in 
the literature, including in ^ [21 El El El El HI El EE EH El and references therein. These variants will be 
discussed in detail in Section 2. In particular, in j7|, the problem of maximizing an objective function 
f (w 1 (■),... ,w d (-)) with / convex and d fixed was considered for combinatorial optimization families in 
general. It was shown, extending earlier results of jSl|H], that if the polyhedra underlying the problem 
have nice edge symmetry (few edge-directions) then the problem can be solved in strongly polynomial 
time. This resulted in polynomial time algorithms for maximizing a convex / for various problems 
including vector partitioning, matroids, and transportation problems with fixed number of suppliers. 
However, the Birkhoff polytope which underlies our bipartite matching problem has exponentially 
many edge-directions (Proposition 12.31 Section 2) and hence the methods of El E] fail- 
Nonlinear bipartite matching is generally intractable, since already for fixed d = 1 (single weight 
function) , the problem of minimizing a family of very simple convex univariate functions f u : R — ► R 
defined by f u (y) '■= (y—u) 2 with u an integer parameter is NP-hard (Proposition ^, 2| part 1, Section 2). 
Therefore, for the most part, the complexity of our results will depend on the unary size of the 
weights, that is, on max In particular, our algorithms will have polynomial complexity for binary 
weights, that is, with wfj S {0, 1} for all i, j, k. In this case, letting Ek be the support of w k for each 
k, the problem becomes that of finding a perfect matching M C N maximizing (or minimizing) 
f(\M D -Ei|, . . . , \M H -Erfl). The problem with binary weights is not easy either: the complexity with 
/ an arbitrary function is unknown for any fixed d > 2; and for variable d it is again NP-hard for 
minimizing the convex multivariate extension of f u above, i.e. the family of functions f u : R rf — > R 
defined by f u (y) '■= Ylt=i(yk ~ u k) 2 and parameterized by u 6 Z d (Proposition 12 .21 part 2, Section 2). 

Clearly, the complexity of the problem depends also on the presentation of the function /: we will 
mostly assume that / is presented by a comparison oracle that, queried on y,z £ Z d , asserts whether 
f(y) < f( z )- This is a broad presentation that reveals little information on the function, making the 
problem harder to solve. In particular, if d is variable, then already for binary weights and maximizing 
a convex /, an exponential number of oracle queries is needed (Proposition 12.21 part 3, Section 2). 

In spite of these difficulties, we are able to provide the following efficient algorithms for the problem: 
in the statements below, oracle-time refers to the running time plus the number of oracle queries. 

Our first theorem provides an efficient algorithm for maximizing convex functions. 

Theorem 1.1 For any fixed d, there is an algorithm that, given any positive integer n, any integer 
weights w 1 ,. . . ,w d , and any convex function f : M. d — > R presented by comparison oracle, solves the 
maximum nonlinear bipartite matching problem in oracle-time which is polynomial in n and max 

A second theorem provides an efficient randomized algorithm for any function. 
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Theorem 1.2 For any fixed d, there is a randomized algorithm that, given any positive integer n, any 
integer weights w 1 ,. . . ,w d , and any function f : R d — > R presented by comparison oracle, solves the 
nonlinear bipartite matching problem in oracle-time which is polynomial in n and max 

We also consider the minimum and maximum nonlinear bipartite matching problems where the 
function / is the l p norm || • \\ p : R d — ► R given by ||y|| p = {Ylt=i \Vk\ p )^ for 1 < p < oo and 
1 1 y | |oo = max^ =1 \yk\- For l p norm minimization we give an algorithm which is polynomial in n and 
max \ uiij\ and determines a (^-approximative solution for any p and a more accurate, \/d-approximative 
solution, for the case of the Euclidean norm p = 2 (Theorem 14. For l p norm maximization we give 
an algorithm which is polynomial even in the bit size of the weights wfj and even if d is variable, and 
determines a dp -approximative solution for any p (Theorem 14.2(1 . 

The article proceeds as follows. In Section 2 we discuss various variants and relatives of the problem, 
survey what is known in the literature about their complexity, and demonstrate the intractability of 
the problem under various conditions. In Section 3 we discuss convex maximization and prove Theorem 
11.11 In Section 4 we discuss approximative norm minimization and maximization and prove Theorems 
14.11 and 14.21 Finally, in Section 5 we discuss randomized optimization and prove Theorem 11.21 

2 Variants and Intractability 

We now discuss various variants and relatives of nonlinear bipartite matching, survey what is known 
(and unknown) about their complexity, and demonstrate its intractability under various conditions. 

First, consider the following related decision problem, asking for the existence of a perfect matching 
attaining specified values under each of the given weight (linear objective) functions w 1 , . . . ,w d . 

Specified multiobjective bipartite matching. Given d,n, weight functions w 1 ,. . . ,w d : N — ► Z, 

and integers u\, ■ ■ ■ , Ud, decide if there is a perfect matching M C N satisfying w k {M) = Uk for all k. 

Chandrasekaran et. al. considered the problem with a single objective w (fixed d = 1) and have shown 
that already this special case is NP-complete This raises the question about the complexity in 
terms of the unary size max|io* -| of the weights. Indeed, even the case of binary weights w\ ^ G {0, 1} 
is not yet understood: for d = 1 it was posed as intriguing and mysterious by Papadimitriou and 
Yanakakis [01 E3) an d the solutions obtained consequently (first by Karzanov 4^ and recently in |12j ) 
are rather sophisticated; for d = 2, the complexity is long open; and for variable d it is NP-complete. 

The following proposition summarizes the known intractability facts about the specified multiob- 
jective bipartite matching problem. We include the short proof for completeness of the exposition. 
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Proposition 2.1 The specified multiobjective bipartite matching problem is NP-complete already un- 
der the following restrictions: (1) fixed d = 1 (single weight function) ; (2) binary weights wfj G {0, 1}. 

Proof. (1) is from |T[ by reduction from subset sum: given integers ao, ai, . . . , a m , define d := 1, n := 2m, 
single weight u> G % nxn by Wij := a, for 1 < i, j < m and Wij := otherwise, and u := ao- Then there 
is a perfect matching M with w(M) = u if and only if there is an / C {1, . . . , m} with ao = Yliei ai ^ 
(2) is by reduction from 3-dimensional matching: given binary n x n x n array x, put d := n, wfj := 
Xjj^, and Uk '■= 1 for all k. Then there is a perfect matching M with w k {M) = u*. = 1 for all if 
and only if there is a binary array y < x with Y^ij Ui,j,k = Yli,k ViJ,k = J2j,k Vi,j,k = 1 for all i,j, k. □ 

A further specialization of the case of binary weig hts wfj E {0, 1} arises when the w k have pairwise 
disjoint supports. This can be formulated as the following particularly appealing "colorful" problem. 

Colorful bipartite matching. Given any bipartite graph G with d-colored edge set E = (+jf =1 Ef, and 
U\, . . . ,Ud, decide if there is a perfect matching MCE 1 containing Uk edges of color k for k = 1, . . . , d. 

This problem is a special case of specified multiobjective bipartite matching with binary weights. To 
see this, note that we may assume G has the same number n of vertices on each side, making it a 
subgraph of K n>n with EON, and Ylt=i u k = n > e ^ se G nas no colorful perfect matching; now, letting 
w k G {0, l} nxn be the indicator of E k for all k, we have that M C N is a perfect matching of K n n 
with w k (M) = Uk for all k if and only if M is a perfect matching of G with \M fl E^\ = Uk for all k. 

For d = 2 (two colors), this problem is sometimes referred to in the literature as the exact matching 
problem: for G = K n ^ n it is polynomial time decidable |UE3; fo r arbitrary bipartite graph G there is 
a randomized algorithm but the deterministic complexity is a longstanding open problem. 

Returning to nonlinear bipartite matching, the next proposition describes its intractability under 
various conditions. By saying that an optimization problem (rather than a decision problem) is NP- 
hard, we mean, as usual, that there can be no polynomial time algorithm for solving it unless P=NP. 

Proposition 2.2 The following hold for the nonlinear bipartite matching problem with data d, n, 
weights w 1 , . . . , w d G Z nxn , and function f : M. d — ► R presented explicitly or by a comparison oracle: 

1. For fixed d = 1 (single weight function) and minimizing the simple convex function f u : R — ► R 
defined by f u (y) '■= (y — u) 2 with u an integer input parameter, the problem is already NP-hard. 

2. For binary weights w\,- G {0, 1} and minimizing the convex function f u : M. d — ► R defined by 
fu(y) '■= Ylt=i(yk~ u k) 2 with u = (ui, . . . ,U(i) an integer vector, the problem is already NP-hard. 
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3. For binary weights wfj and maximizing a convex f presented by comparison oracle, exponentially 
many oracle queries are needed and hence the problem is not solvable in polynomial oracle-time. 

Proof. 

1. Given weight w and integer u, there is a perfect matching M with w{M) = u if and only if the 
minimum value f u (w(M)) = (w(M) — u) 2 of a perfect matching M under f u is 0. So even computing 
the optimal objective function value enables deciding the NP-complete problem (1) of Proposition 12. 11 

2. Analogously to the proof of part 1 above: given binary weights w l , . . . ,w d and vector u = (u\, . . . , Ud), 
there is a perfect matching M with w k (M) = Uf- for all k if and only if the minimum objective value 
f u (w 1 (M), . . .,w d (M)) = £fc =1 (w fc (M) - u k ) 2 of a perfect matching M under f u is 0. bo even com- 
puting the optimal objective value enables deciding the NP-complete problem (2) of Proposition 12.11 

3. Let d = n 2 , define binary weights w r,s G Z" xn for 1 < r, s < n, with iup? := 1 if (i,j) = 
(r,s) and = otherwise, and let / : M. d = M nxn — ► R be any function. Then for any matrix 
x £ M nxn we have f(w 1 ' 1 x, . . . , w n,n x) = f(xi t i, . . . , x njU ) = f(x). Since the permutation matrices 
(which correspond to perfect matchings) are convexly independent, any assignment of values to the n! 
permutation matrices can be extended to a convex function / on R nx ". Thus, to find the permutation 
matrix maximizing /, the oracle presenting / must be queried on all n! permutation matrices. □ 

More generally, the nonlinear combinatorial optimization problem is the following: given positive 
integers d, n, a family T of subsets of a ground set {1, . . . , n}, integer weights w 1 , . . . , w d on {1, . . . , n}, 
and an arbitrary function / : M rf — > K, find F G T maximizing (or minimizing) f(w 1 (F), . . . , w d (F)). 
In the maximization problem with / convex and d fixed was studied. It was shown that if the 
number of edge-directions of the polytopes P T := convjli? : F G J 7 } (where lp G {0,1}™ denotes 
the indicator of F) is polynomial in n for a class of families presented by membership oracles then the 
problem over families T in that class can be solved in strongly polynomial oracle-time. This unified 
and extended earlier results of [HI [HI and yielded polynomial time algorithms for convex maximization 
for various problems including vector partitioning, matroids, and transportation problems with fixed 
number of suppliers. However, for bipartite matching, which is the combinatorial optimization problem 
over the family T C 2^ of perfect matchings in K n>n , the underlying polytope is the Birkhoff polytope 

= u n := {ie»f : ^2xij = l, ^x id = 1} 

i j 

which, as we next show, has exponentially many edge-directions, making the methods of [HI 13 El fail- 
Proposition 2.3 The Birkhoff polytope Tl n has precisely ^Ylk=2 (T) 2 ^^—^- — nCa) edge- directions. 
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Proof. Every edge-direction of IT™ is a nonzero minimal-support matrix x G W nxn with zero row-sums 
and column-sums, and hence (up to scalar multiplication) is the matrix xc of some circuit C C N 
of K n>n , having values ±1 alternating along the edges of the circuit and elsewhere (see e.g. [7.). 
We claim that each such circuit matrix xc is an edge-direction. To see this, let C = C + bj C~ be 
the partition of alternating edges of C and let D be a matching in K n ^ n which perfectly matches all 
vertices not in C . Let x + and x~ be the permutation matrices which are the indicators of the perfect 
matchings M + := C + U D and M~ := C~ U D of K n ^ n . Define a binary weight matrix w as the 
indicator of CUD. Then wx + = wx~ = n whereas wx < n for any other permutation matrix x. Thus, 
wx attains its maximum over IP precisely at the two vertices x + and x~ and hence [x + , x~] is an edge 
and the difference xc = x + — x~ is an edge-direction. Now, for each k > 2, the number of 2/c-circuits 
of K n ^ n is known and easily seen to be \ {^) 2 k\{k — 1)! (see [7j) and hence the Proposition follows. □ 

Proposition ^ . 31 shows that, while the methods of [HJEllHj do apply for transportation problems with 
fixed number of suppliers, they fail for bipartite matching which is the simplest possible transportation 
problem - albeit, with variable numbers of suppliers and consumers - and do not lead to a polynomial 
time algorithm even for maximizing a convex / with d fixed. This state of things, along with the easy 
solvability of the standard linear bipartite matching problem, make the nonlinear problem for bipartite 
matching particulary intriguing, and is part of our motivation in raising and studying it herein. 

3 Deterministic Convex Maximization 

In this section we discuss the maximum nonlinear bipartite matching problem for convex functions 
/ : R d — ► M presented by comparison oracles. We start with some definitions. Here we will be 
working with matrices rather than graphs and matchings, so the weights are now integer matrices 
w 1 , . . . ,w d G Z nxn , and the solutions are permutation matrices, which are well known to be precisely 
the vertices of the Birkhoff polytope of bistochastic matrices (with M + the nonnegative reals), 

IT = {x G R n + * n . J2 XiJ = 1, Y^Xij = 1} . 

Given weights w 1 , . . . , w d , define a projection w : R nxn — > M. d mapping matrices x to vectors w ■ x, 
w ■ x := (w l x, . . . , w d x) = ( y~] wjjXij, . . . , wfjXij) . 

i,3 h3 

Define the multiobjective polytope (corresponding to w 1 , . . . , w d ) to be the projection of 11™ under w, 
n™ := {wx = (w 1 x,...,w d x) : x G IT} C M. d . 
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Finally, define the fiber of any point y £ M. d to be the polytope H n Clw~ 1 (y) consisting of those matrices 
in the Birkhoff polytope that are projected by w onto y. Thus, a point y is in II™ if and only if its fiber 
is nonempty; the following lemma asserts that these equivalent conditions can be decided efficiently 

Lemma 3.1 There is a polynomial time algorithm that, given d,n,w l , . . . ,w d , and integer y £ Z d ; 
either asserts y £" II™ and II™ (~l w~ 1 (y) = or asserts y £ II™ and returns a vertex x ofH n n w~ 1 (y). 

Proof. The fiber of any y = (yi, . . . , yd) is the polytope given by the following inequality description, 

ffnr'fe) = {££11™ : (w 1 x,...,w d x) = (y 1 ,...,y d )} 

= {:r£lR™ x ™: ^^ = 1, ^^ = 1, w k x = y k } , 

i j 

so linear programming allows to efficiently compute a vertex of the fiber or assert that it is empty. □ 
This lemma implies in turn that the multiobjective polytope II™ can be constructed efficiently. 

Lemma 3.2 For any fixed d, there is an algorithm that, given n and w 1 , . . . , w d £ Z™ x ™ ; computes the 
vertex set vertili 1 ^) of the multiobjective polytope II™ in time which is polynomial in n and max \wf A. 

Proof. Let u : = max|ti;^|. Then for any permutation matrix x and its projection y = w ■ x we 
have \yk\ = \w k x\ < nu and therefore y lies in the grid {0, ±1,... , ±nu} d . Since each vertex y of 
II™ is the projection y = w ■ x of some vertex x of II™, which is a permutation matrix, we have 
vert(II™ ) C {0, ±1, . . . , ±nu} d . For each of the (2nu + l) d grid points y £ {0, ±1, . . . , ±nu} d , apply 
the algorithm of Lemma 13.11 to check if y £ II™ , and obtain Y := {0, ±1, . . . , ±nu} d n II™ . 

We then have that vert (II™ ) C Y C II™ and therefore the multiobjective polytope II™ is the convex 
hull of Y. Since convex hulls can be computed in polynomial time for any fixed dimension d, we can 
efficiently construct II™, that is, determine all its vertices (and more generally all its faces). □ 

Lemma f3 . 1 1 shows that for any y £ Z d it is possible to check efficiently if y is the projection y = w-x 
of some bistochastic matrix x £ II™, and to find such an x if one exists. We will need also to consider 
the integer analog of this problem: given y £ Z d , is y the projection y = w ■ x of some permutation 
matrix x £ vert(II™), and if it is, can we find one such x efficiently? but this problem is precisely the 
specified multiobjective bipartite matching problem: there is a permutation matrix x with y = w ■ x 
if and only if there is a perfect matching M with w k (M) = y k for k = 1, ... ,d. Unfortunately, as 
explained in Section 2, the complexity of this problem is open even for fixed d = 2. The difficulty is 
that the fiber II™ n w^ 1 (y) of y is not necessarily an integer polytope and it may have some fractional 
(bistochastic) matrices and some integer (permutation) matrices x as its vertices. 

Fortunately, as the next lemma shows, the fibers of vertices of II™ are better behaved. 
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Lemma 3.3 Let y be any vertex of the multiobjective polytope II™ . Then the fiber II™ n w^ 1 (y) of y is 
a nonempty integer polytope all of whose vertices are permutation matrices. Thus, the polynomial time 
algorithm of Lemma \3.1\ applied to y G vertijl 7 ^) returns a permutation matrix x satisfying y = w ■ x. 

Proof. It is well known and easy to see that if Q is the image of a polytope P under an affine map a, 
then the preimage P n a^ 1 (F) = {x G P : a(x) G F} of any face F of Q is a face of P. Thus, if y is 
a vertex of II™ then its fiber II n n u? _1 (y), which is the preimage under the map w of the face {y} of 
II™ , is a face of Il n . Therefore, the vertices of the nonempty fiber of y, one of which will be returned 
by the algorithm of Lemma 13.11 are precisely the vertices of IP which are contained in that fiber. □ 

We can now prove our first theorem, providing an efficient algorithm for convex maximization. 

Theorem 11.11 For any fixed d, there is an algorithm that, given any positive integer n, any integer 
weights w 1 , . . . , w d , and any convex function / : M. — > IR presented by comparison oracle, solves the 
maximum nonlinear bipartite matching problem in oracle-time which is polynomial in n and max • \ . 

Proof. Since / is convex on M d and f(w 1 (-),...,w d (-)) is convex on R™ x ™, and the maximum of a 
convex function over a polytope is attained at a vertex of the polytope, we have the following equality, 

max {f(w 1 x, . . . , w d x) : x £ vert (If 1 )} = max{f(w 1 x,...,w d x): x G Il n } 

= max{/(y): yen™} = max{/(y) : y G vert (II™)} . 

Apply the algorithm of Lemma 13.21 and compute vert (II™). By repeatedly querying the comparison 
oracle of /, identify a vertex y* G vert (II™) attaining maximum value f{y). Now apply the algorithm 
of Lemma l3.1l to y* and, as guaranteed by Lemma l3,3| obtain a permutation matrix x* in the fiber 
of y*, so that y* = w ■ x* = {w 1 x* , . . . ,w d x*) and f(w 1 x*, . . . ,w d x*) = f(y*). Since y* attains the 
maximum on the right-hand side of the equation above, x* attains the maximum on the left-hand 
side. Thus, the perfect matching of K n ^ n corresponding to the permutation matrix x* is optimal. □ 

The most time consuming part of the algorithm underlying Theorem 11.11 is the repeated use of 
linear programming for testing fibers of points in the grid {0, ±1, . . . , ±nu} d to construct vert(II™ ). 
There are various ways of improving the algorithm in practice, but they do not seem to improve the 
worst case complexity. We now describe such a variant of the algorithm which will usually be much 
faster since it will typically test the fibers of some but not all points in the grid. 

A variant of the convex maximization algorithm. 

1. Find the smallest grid containing vert (II™ ) by solving, for k = 1, . . . , d, the two linear programs 
Sfc := mm{w k x : Xij = Xjj = 1, x > 0} , ifc := max{w k x : Xij = Xjj = 1, x > 0} ; 

i j i j 



8 



then vert (II™ ) is contained in the grid Z := {y G Z d : < y^ < t^, k = 1, . . . , d}. 

2. By repeatedly querying the comparison oracle of /, order the grid points by nonincreasing value 
under / and label them y 1 , . . . , y\ z \ so that Z = {y 1 , . . . , and f{y l ) > • • • > 

3. Apply the algorithm of Lemma 13.11 to test the fiber of each y l in order, until the first y k for 
which the vertex x* of its fiber W 1 f\w^ 1 {y k ) returned by the algorithm is a permutation matrix. 

4. Output the perfect matching of K n ^ n corresponding to the permutation matrix x* . 

We claim that x* is an optimal solution to the maximum convex bipartite matching problem. Indeed, 
note that /* := max{/(y) : y G II™ } = max{/(y) : y G vert(Il™ )} equals the optimal objective function 
value max{/(w 1 x, . . . ,w d x) : x G vert(n n )} (see proof of Theorem 1 1.1 (I ; let y m G vert(n^) C Z be a 
vertex achieving that maximum value f(y m ) = /*; by Lemma l3,31 the algorithm of Lemma l3.1l applied 
to y m G vert(n^) returns a permutation matrix and so k < m; this implies /* > f{y k ) > f{y m ) = f* 
and hence /* = f{y k ) = f(w 1 x*, . . . , w d x*); therefore x* achieves the optimal objective function value. 



We end this section with an example of a maximum convex bipartite matching problem, demon- 
strating all notions and algorithms discussed above, some of which will be also used in later sections. 



Example 3.4 Consider the maximum convex bipartite matching problem with the following data: 
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By solving the linear programs minimizing and maximizing w k x over II 4 for k = 1,2 (step 1 of the 
algorithm above) we get si = S2 = 0, t\ = 3, t% = 4, and so the smallest grid containing vert(n 4 ) is 
Z := {y G Z 2 : < y 1 < 3, < y 2 < 4} C {0, ±1, . . . , ±4} 2 which contains 20 points. Figure □ below 
depicts this grid and indicates the objective function value f(y) = yf + y 2 of each grid point. Ordering 
the points by decreasing value under / (step 2 above) we get y 1 = (3, 4) , y 2 = (2, 4) , ... , y 20 = (0, 0). 
Testing fibers of the y l in order (step 3 above), the fiber of y 1 is found empty whereas the fiber of y 2 
is nonempty and is the first for which the algorithm of Lemma 13.11 returns a permutation matrix 
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Thus, the corresponding matching M* := {(1, 1), (2, 4), (3, 3), (4, 2)} of K4 4 is an optimal solution. 

Figure H also shows the multiobjective polytope 11^ and its vertex set vert(II^) computed by the 
algorithm of Lemma l3.2l blue circles are non- vertex grid points in 11^ and green diamonds are vertices 
of 11^,. The optimal point y 2 = (2, 4) which is found either by the algorithm above or by the algorithm 
of Theorem ll.ll is the vertex of 11^ attaining maximum value under / and is a red square. Of particular 
interest is the blue point y = (1, 2) whose fiber 11" C\w~ 1 (y) is a non-integer polytope with 30 vertices 
(more than the 24 of the Birkhoff polytope upstairs!), all of which are fractional, such as 
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indicating the difficulty of the specified multiobjective and colorful bipartite matching problems. 
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Figure 1: Example 13.41 and the polytope II^ 



4 Approximative Norm Optimization 

Consider any discrete optimization problem with a finite set S of feasible solutions and nonnegative 
objective function g : S — ► M+ to be minimized or maximized and let s* S S be any optimal solution. 
Then an r -approximative solution is any feasible solution s £ S satisfying ^g(s*) < g(s) < rg(s*). 



10 



In this section, building on the tools and results of Section 3, we provide approximative algorithms 
for the minimum and maximum nonlinear bipartite matching problems where the function / is the 
l p norm || • || p : R d — > E given by \\y\\ p = (J2k=i \Vk\ p )* for 1 < p < oo and \\y\\oo = maxf =1 \y k \. 
To keep our results general and allow treatment of fractional and even nonrational p, we will still 
assume that / is presented by a comparison oracle. Of course, for the most common values p = 1, 2, oo 
such an oracle is realizable in time polynomial in the rest of the data; moreover, for any integer p, by 
computing and comparing the integer valued p-th. power ||y||p of the norm instead of the norm itself, 
such an oracle is realizable in time polynomial in the rest of the data and [logp]. 



4.1 Minimization 



The following theorem provides an efficient approximative algorithm for minimizing the L norm. 



Theorem 4.1 For any fixed d, there is an algorithm that, given any n, any 1 < p < oo, and any non- 
negative integer weights w 1 , . . . ,w d , determines a d- approximative solution to the minimum nonlinear 



,, in oracle-time which is polynomial in n and maxwj 



bipartite matching problem with / = || 
For p = 2 (Euclidean norm), the algorithm determines a more accurate, Vd -approximative, solution. 



Proof. The algorithm is the following: apply the algorithm of Lemma 13.21 and construct the ver- 
tex set vert (II™) of the multiobjective polytope. Using the comparison oracle of / identify a vertex 
y 6 vert (II™) attaining minimum value ||y|| p . Now apply the algorithm of Lemma 13.11 to y and, as 
guaranteed by Lemma 13.31 obtain a permutation matrix x in the fiber of y, so that y = w ■ x. Output 
the perfect matching of K nn corresponding to the permutation matrix x. 

We now show that this provides the claimed approximation. Let x* be the permutation matrix 
corresponding to an optimal perfect matching and let y* := w-x* be its projection. Let y' be a point on 
the boundary of II™ satisfying y' < y* . By Caratheodory's theorem (on the boundary) y' is a convex 
combination y' = £[ =1 ^iV z °^ some r < d vertices of II™ and hence Af = max Aj > - > -j for some t. 

Since the weights w k are nonnegative we find that so are y' and the y l and hence we obtain 



f(w l x, . . . , w d x) 



\y\\p 



< 



\y% < d\ t • |||/* 



d- || V| 



d ■ 



Ev 



i=l 



d-\\y'\\ P < d-\\y* 



d- f(w 1 x*,...,w d x*) 



This proves that x provides a d- approximative solution for any 1 < p < oo. Now consider the case of 
Euclidean normp = 2. By Cauchy-Schwartz, 1 = (£T =1 1-A;) 2 < YJi=i ! 2 E[=i >H = r £ A? < d£ A?. 
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Find s with \\y s \\ p = min \\y l \\ p and recall that the y l are nonnegative. We then have the inequality 

r r 

f\w l x,...,w d x) = \\y\\ 2 2 < \\y s \\ 2 2 < (d£ A?) ■ \\y s \\l < dJ^^lWvll 

i=i i=i 

= d- Us/Hi < d-||i/*||l = d-/VV,...,«;V) 



i=i 



2 



which proves that in this case, as claimed, x provides moreover a v^-approximative solution. □ 



4.2 Maximization 

The following theorem provides an approximative algorithm for maximizing the l p norm, that runs in 
time which is polynomial even in the bit size of the weights wfj and even if d is variable. 

Theorem 4.2 There is an algorithm that, given any d, any n, any 1 < p < oo, and any nonnegative 
integer weights w , . . . , w , determines a dp -approximative solution to the maximum nonlinear bipartite 



ma 



tching problem with / = || ■ || p , in oracle-time which is polynomial in d ,n, and maxflogu;^ ■] . 



Proof. The algorithm is the following: for k = 1, . . . , d solve the linear programming problem 

max{u; fc x : Xjj = 1, Xjj = 1, x > 0} , 

» i 

obtain an optimal vertex x k of IT 1 , and let y k := w ■ x k be its projection. Using the comparison oracle 
of / find r with \\y r \\ p = max^ =1 ||y fc || p . Output the perfect matching of K n ^ n corresponding to x r . 

We now show that this provides the claimed approximation. Let s satisfy ||y s ||oo = max^ =1 ||y fc ||oo- 
First, we claim that any y £ IT™ satisfies ||y||oo < Halloo- To see this, choose any point x G H n r\w~ 1 (y) 
in the fiber of y so that y := w ■ x, let t satisfy yt = \\y\\oo = max^ =1 y^, and recall that the w k and 
hence the y k are all nonnegative. Then, as claimed, we get 



Mloo = yt = w l x < maxjai'i : x £ Il n } = w t x t = y\ < ||y*||oo < ||y s ||oo 



Let x* be an optimal permutation matrix and let y* := w ■ x* be its projection. Consider first the 
case p = oo. Then ||y r ||oo = max^ =1 ||y fc ||oo = ||?/ s ||oo and hence, by the claim just proved, we have 

_/Vx*,...,™V) = Hyloo < ||2/ s ||oo = Halloo = /(™V,...,™V) < /(wV,..,^*) . 
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Therefore equality holds all along and x r provides an exact optimal solution, or in other words, a 
1-approximative solution, agreeing with the statement of the theorem with d~ = 1 for p = oo. Next, 
consider the case of any 1 < p < oo. Then we have the following inequality which completes the proof, 

d 

nw 1 x\...,w d x*) = \\y*\\v = El^l" ^ d-\\yToo < d-\\yTc 

k=i 

d 

< d-J2\v S k\ P = d-\\yX < d-\\y r \\l = d-f{w 1 x s ,...,w d x s ).n 
fc=i 

5 Randomized Nonlinear Optimization 

In this section we provide a randomized algorithm for nonlinear bipartite matching for any function 
/ : M. d — ► R presented by a comparison oracle. By this we mean an algorithm that has access to a 
random bit generator, and on any input outputs the optimal solution with probability at least half. 

By adding to each wfj a suitable positive integer v and replacing the function / by the function 
that maps each y G M. d to f(y\ — nv, . . . , — nv) if necessary, we may and will assume without loss 
of generality throughout this section that the given weights are nonnegative, w 1 , . . . ,w d € N nxn . 

Recall that vert(Il ra ) is the set of n x n permutation matrices and let Y := {w ■ x : x £ vert(II n )} 
be the set of all projections y = w • x = (w 1 x, . . . , w d x) G N d of permutation matrices x. 

In this section we will be working with polynomials with integer coefficients in the n 2 + d variables 
Ojj, i, j = 1, . . . , n and b^, k = 1, . . . , d. Define an n x n matrix A whose entries are monomials by 

d k 

A i,3 : = a i,jT[ b 7' j ' i,j = l,...,n . 

k=l 

For each matrix x £ f$ nxn and vector y G N n , the corresponding monomial is 

n n d 

<rv := 1I1K; ll'f • 

1=1 j=l k=l 

For each permutation matrix x let sign(x) = ± denote the sign of the corresponding permutation. 
Finally, for each y 6 N d define the following polynomial in the variables a = (a^j) only, by 

g y (a) := ^ {sign(x)a x : x £ vert(II n ), w • x = y} 

We then have the following identity expanding the determinant of A in terms of the g y (a) , 

det(A) = Yl sign(x)n<r = E ^(^ X b WX = E 9y(a)b' y . 

a;evert(n") i,j xevert(n™) j/Gy=w)-vert(n™) 
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Next we consider integer substitutions to the variables aij. Under such substitutions, each g y (a) 
becomes an integer and det(A) = J2 y eY 9y( a )b v becomes a polynomial in the variables b = only. 
Given such a substitution, let Y := {y £ Y : g y (a) ^ 0} be the support of det(^4), that is, the set of 
exponents of monomial b y appearing with nonzero coefficient in det(A). 

The next proposition concerns substitutions of independent identical random variables uniformly 
distributed on the set of integers {1, 2, . . . , s}, under which Y becomes a random subset of Y. 

Proposition 5.1 Suppose that independent identical random variables uniformly distributed on the 
set {1, 2, . . . , s} are substituted for the a^j and letY = {y£Y : g y {a) ^ 0} be the random support of 
det(.A). Then, for every y &Y = {w ■ x : x £ vert(n. n )} , the probability that y ^Y is at most —. 

Proof. Consider any y £ Y and consider g y (a) as a polynomial in the variables a = (ajj). Since y = w-x 
for some permutation matrix, there is as least one term sign(x)a x in g y {a). Since distinct permutation 
matrices x give distinct monomials a x , no cancellations occur among the terms sign(x)a x in g y (a). 
Thus, g y (a) is a nonzero polynomial of degree n. The claim now follows from a lemma of Schwartz 
jllj stating that the substitution of independent identical random variables uniformly distributed on 
{1, 2, . . . , s} into a nonzero multivariate polynomial of degree n is zero with probability at most j. □ 

The next lemma shows that, given a,j, the support Y of det(yl) is polynomial time computable. 

Lemma 5.2 For any fixed d, there is an algorithm that, given n, w 1 , . . . ,w d £ pp xn ; and substitutions 
aij £ {1, 2, . . . , s}, computes Y = {y £ Y : g y (a) ^ 0} in time polynomial in n, maxw^j and [log s] . 

Proof. For each y, let g y := g y (a) be the fixed integer obtained by substituting the given integers a^j. 
Put u = n ■ rnaxu),^ and Z = {0, 1, . . . ,u} d . Then Y C Y C Z and hence det(yl) = J2 y ez 9y^ V ls a 
polynomial in d variables b = (bk) involving at most \Z\ = (u+l) d monomials. For t = 1,2,..., (u+l) d 
consider the substitution b^ := 1 ; k = l,...,d. Let A(t) be the integer matrix obtained from 

A by this substitution along with the substitution of the given Gtjj. Then each entry of A(t) satisfies 

m<j = ^n^" 1 )^ < s n(«« +i ) rf ) (tt+ri ) f ^ S (u+i) d ^ d+i 

and hence its bits size 1 + log A(t)ij = 0(u d+1 log(su)) is polynomially bounded in n, maxwfj, [log s] . 
Therefore the integer number det(A(t)) can be computed in polynomial time by Gaussian elimination. 
So we obtain the following system of (u + l) d equations in (u + l) d variables g y , y £ Z = {0, 1, . . . , u} d , 

det(A(t)) = Y.9v\l h l k = J2 t ' ELlVkiU+1)k " ■ 9v . t = 1,2, . . . , (u + l) d . 

yez k=l yez 
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As y = (yx, . . . , yd) runs through Z, the sum Ylt=i Vk{ u + ^) k ~ 1 attains precisely all \Z\ = (u + l) d 
distinct values 0, 1, . . . , (u + l) d — 1. This implies that, under the total order of the points y in Z by 
increasing value of Ylt=i Dk{u + ^) k ~ l , the vector of coefficients of the g y in the equation corresponding 
to t is precisely the point (t ^ 1 , . . . , t(" +1 ) d ~ 1 ) on the moment curve in M z ~ R( u+1 ) d . Therefore, the 
equations are linearly independent and hence the system can be solved for the g y = g y (a) and the 
desired support Y = {y £ Y : g y (a) ^ 0} of det(^4) can indeed be computed in polynomial time. □ 

We are now in position to prove Theorem ll.2l By a randomized algorithm that solves the nonlinear 
bipartite matching problem we mean an algorithm that has access to a random bit generator and on 
any input to the problem outputs a perfect matching which is optimal with probability at least a half. 
The running time of the algorithm includes a count of the number of random bits used. Note that 
by repeatedly applying such an algorithm several times and picking the best perfect matching, the 
probability of failure can be decreased at will; in particular, repeating it n times decreases the failure 
probability to as negligible a fraction as while increasing the running time by a linear factor only. 

Theorem 11.21 For any fixed d, there is a randomized algorithm that, given any positive integer n, 
any integer weights w 1 , . . . , w , and any function / : M. d — ► R presented by comparison oracle, solves 
the nonlinear bipartite matching problem in oracle-time which is polynomial in n and max - |. 

Proof. As explained in the beginning of this section, we may and will assume that the w k are non- 
negative. First we claim that, with probability at least 1 — Jj-, we can compute the optimal objective 
function value of the nonlinear bipartite matching problem. To see this, note that the optimal value 
equals max{/(y) : y £ Y} where Y = {w ■ x : x £ vert(n ra )} as before, and let y* G Y be a 
point attaining f(y*) = max{/(y) : y G Y}. Now, using polynomially many random bits, draw in- 
dependently and uniformly distributed integers from {1,2, . . . ,2n 2 } and substitute them for the ajj. 
Next compute Y = {y £ Y : g y (a) ^ 0} using the algorithm underlying Lemma 15.21 and determine 
max{/(y) : y £ Y}. By Proposition 15.11 with probability at least 1 — we have y* £ Y in which 
event max{/(y) : y £ Y} = max{/(y) : y £ Y} is indeed the optimal objective function value. 

Next, suppose that M C N = ■ 1 < i,j < n} is any (not necessarily perfect) matching of 

K n>n . Then we can also compute, with probability at least 1 — the maximum objective function 
value among perfect matchings MUL containing M. To see this, let m := n — \M\ and consider 
the subgraph G of K n-n induced by the vertices not matched under M. Then G is isomorphic to 
Km,m and we have a naturally induced nonlinear bipartite matching problem on G, where the new 
weight functions w k are simply the restrictions of the w k to the edges of G, and the new functional 
/ on M. d is defined by f(y±,..., y^) '■= f(yi + w 1 (M), . . . , y^ + w d (M)). Then the objective function 
value f(w 1 (M U L), . . . , w d (M UL)) of any perfect matching M U L of K n ^ n in the original problem 
equals the objective function value f(w l (L), . . . , w d (L)) of the perfect matching L of G in the induced 
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problem. Since max wfj < maxwfj and m < n we can compute, with probability at least 1 — in 
time polynomial in n and max \ w^j\, the optimal objective function value of a perfect matching of G 
by the algorithm of the paragraph above applied to G, where the randomized substitutions are taken 
from {1, 2, . . . , 2mn} (and not from {1, 2, . . . , 2m 2 } which would give smaller probability of success). 
This value is the maximum objective function value among perfect matchings M U L containing M. 

We claim that the following procedure constructs a perfect matching M of K n>n which is optimal 
with probability at least \. Start with M := and i := 1. While i < n iterate: for each edge such 
that j is not matched under M, use the algorithm of the previous paragraph to obtain the maximum 
objective function value of a perfect matching of K n>n containing M U {(i, j)}; let ji be the smallest j 
for which this value is maximal; update M := M L) increment i and repeat. Output M. 

To prove the claim, let M* = {(1, n), (2, r^), . . . , (n, r n )} be the lexicographically first optimal per- 
fect matching of K nn , that is, the one such that for any other optimal M 1 = {(1, s±), (2, s%), . . . , (n, s n )} 
there is an index 1 < h < n such that r» = Sj for all i < /i and < s^. For i = l,...,n let 
Ei be the random event that after the completion of iteration i of the above procedure we have 
M = {(l,ri), . . . , (i,ri)}. We prove by induction on i that Pr(E'j) > (1 — ^)*. For the basis note 
that Ei is the event that the randomized algorithm used during the first iteration computes correctly 
the maximum objective function value of a perfect matching containing {(l,ri)}, having probabil- 
ity at least 1 — For the inductive step note that Pr(.Ej+i|.Ej) is the probability that, given that 
after iteration i we have M = {(1, r\ ),..., (i, n)}, the randomized algorithm used during iteration 
i + 1 computes correctly the maximum objective function value of a perfect matching containing 
{(1, ri), . . . , (i, rj), (i + 1, rj+i)}, which is again at least 1 — as Ei + \ C Ei, the induction follows by 

Pr(^ +1 ) = Pr(£ m n^) = Pr(^ +1 |^)Pr(^) > (1 - -^)(1 - = (1 - . 

Now, the probability that the perfect matching M output by the procedure above is optimal is no 
smaller than the probability that M equals the lexicographically first optimal perfect matching M* 
which is precisely Pr(.E n ) and hence at least (1 — > ^ as desired. This completes the proof. □ 
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